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Abstract 

A class of loop diagrams in general relativity appears to have a 
behavior which would upset the utility of the energy expansion for 
quantum effects. We show through the study of specific diagrams that 
cancellations occur which restore the expected behaviour of the energy 
expansion. By considering the power counting in a physical gauge we 
show that the apparent bad behavior is a gauge artifact, and that the 
quantum loops enter with a well behaved energy expansion. 
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I Introduction 



Loop calculations in general relativity are readily interpreted using the tech- 
niques of effective field theory fj], H]. As in all effective field theories, the 
utility of such calculations is tied to an expansion in powers of the energy or 
inverse distance. In chiral theories, Weinberg || has provided an important 
theorem which states that diagrams with increasing numbers of loops con- 
tribute to an amplitude with increasing powers of the energy, with each extra 
loop adding an extra factor of E 2 . For example, if one is working to order E A 
accuracy one needs to include only one loop diagrams. While pure gravity 
behaves exactly in the same way, if we try a simple extension of this same 
argument to gravity interacting with matter, we will see in Sec. |I| that the 
desired behavior is not obtained. There is a class of diagrams which appears 
to have Gm 2 as the expansion parameter. This would upset the utility of the 
energy expansion. The purpose of this paper is to explore this problem and 
to see if it obstructs the energy expansion. 

The desired expansion parameters for quantum corrections in an effective 
theory of gravity is Gq 2 ~ ^f, such that at low energies/long distances the 
higher order loop effects are suppressed with respect to tree diagrams and 
low order loops. Thus we can obtain predictions to a given order with a finite 
amount of calculation. General relativity also contains the classical expansion 
parameter Gmq ~ — which represents the nonlinearities of the classical 
theory. This can be found in the loop expansion from the nonanalytic terms 
of the form Gq 2 ^J~^p. However, Gm 2 as an expansion parameter is a major 
problem. In the first place, the mass can be extremely large in units of the 
Planck mass (e. g. m — Afg un ) so that Gm 2 can be a number very much 
larger than unity. In addition if we restore factors of h, this dimensionless 
combination goes like ^P-. The classical limit h ^ would be upset by 
corrections of this form. 

We will see that the apparent difficulty with the loop expansion appears 
to be a gauge artifact. When calculating in harmonic gauge, where the power 
counting is first discussed in Sec. ||, there occur cancellations between indi- 
vidual diagrams, cancelling the bad behavior. We detail the calculation for 
the box and crossed box diagrams. Part of the problem is due to the occur- 
rence of both classical and quantum effects in the same Feynman diagram, 
when treated in covariant gauges. This suggests that separating the classi- 
cal physics from the physical quantum (transverse and traceless) degrees of 
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freedom will improve the power counting. For the interaction of two nearly 
static masses, we show that this is in fact the case. 

The organization of this paper is as follows. In Sec. |TJ we make a naive 
generalization of the Weinberg power counting theorem and isolate those 
diagrams which appear to give a problem. In order to explore this without 
all the tensor indices of gravitons, we introduce a scalar toy model with the 
same behavior in Sect. Ill in order to see how cancellations occur. Sec. |IV| 
applies the lessons so learned to the gravitational interaction. Sec. |V| was 
devoted to development of the power counting scheme in a physical gauge, 
and to the interpretation of the apparent problem as a gauge artifact. We 
end with some concluding comments. 



II Power counting in covariant gauges 





~ Kq 2 ~ n 2 q 2 

(a) (b) 

Figure 1: Three and four graviton couplings. 

We are interested in treating powers of energies and masses in vertices and 
propagators in order to determine the overall energy dependence of a given 
multiloop diagram. The mass of the matter field is not a small parameter, 
but we can treat the external three-momenta as small if we are working at low 
enough energies. Let us review the Feynman rules, and extract the essential 
dependence of the vertices. Starting from the action 

S= [d i x^--^-R (II.l) 

with k 2 = 32tcG, we expand this metric 

9nu = Vvv + (II.2) 
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where h^ v is the fluctuating field. Expanding \R in powers of h^ u we see that 
a term which involves n graviton fields, i. e. h n , carries a coupling constant 
K n ~ 2 . Since the curvature is second order in derivatives, all terms emerging 
from the Einstein action will be of order q 2 . Thus the triple graviton coupling 
of Fig. [I] is of order nq 2 , while the four graviton vertex is of order K 2 q 2 etc. 
The matter fields couple to gravitons through T^, which for a scalar field 
has matrix elements 

< p' | | p >= p^p' u + p v p^ - -g^ (p-p' - m 2 ) (II.3) 

with p^ = (^y/m 2 + p 2 | pj . Treating the mass as a large parameter leads 

to a one graviton vertex (see Fig. which behaves as urn 2 while the two 
graviton diagram is of order K 2 m 2 etc. 



2 2 2 

nm Km 
(a) (b) 

Figure 2: Matter- graviton couplings from T^ v . 

The graviton propagator, like all massless boson propagators, scales as 
^2. The matter field propagator requires a bit more explanation. Because 
we are dealing with an effective theory at low energies, we need not consider 
loops of heavy matter fields. These loops have already been integrated out 
in order to define the low energy effective theory. However, we do need to 
consider matter fields which appear as external states and which propagate 
through a given diagram interacting with each other and with gravitons. The 
explicit form of the propagator is 

Dip + q)= . - ! 2 j = ^ 2 \, 2 (H.4) 

[p + q) — m z 2p ■ q + q z + [p z — m z ) 

where p is the momentum that the matter field has as an external particle, 
and q is the momentum which has been added to it through interactions 
with gravitons (internal or external). The external momentum is on shell 
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(p 2 — m 2 = 0) so that the matter propagator is counted as a factor of 
Note that if we had chosen a different normalization for our matter fields 
(e. g. a nonrelativistic normalization such that T 00 ~ m and D(q) ~ ~ ^ 
both the vertices and propagators change in a way that compensates each 
other, leading to the same counting rules as in our normalization. 

Before giving the general power counting theorem, let us illustrate the idea 
with two specific examples, one of which illustrates the "good" behavior and 
one which shows the problem. First consider graviton-graviton scattering, 
whose overall matrix element is dimensionless. At lowest order we have a 
^2 factor from the coupling in the Einstein action, one of k 4 from the four 
graviton fields and one q 2 because the Einstein action involves two derivatives. 
This leads to an overall matrix element 

Mtree ~ «V- (H.5) 





K 2 q 2 



4 4 

~ k q 



(a) 



(b) 



Figure 3: Sample four-graviton interaction diagrams to one loop, illustrating 
the expected behavior in the energy expansion. 



If we try to iterate this vertex to produce the one loop diagram of Fig. ^b 
we obtain schematically 

Ml °°" ~ J {2nY l 2 (l-q) 2 (IL6) 

where pi,P2,q are various combination of external momenta. If this loop 
integral is regularized dimensionally, which does not introduce powers of 
any new scale, the integral will be represented in terms of the exchanged 
momentum to the appropriate power. Thus we have 

m 1oop ~ (ii.7) 
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where again q represents some combination of external momenta. [There may 

2 

also be logarithms of where \i is the usual scale introduced in dimensional 
regularization.] In this case adding a loop has generated an effect which is 
higher order in the energy expansion. The expansion is in terms of powers 
of K 2 q 2 . 



(a) (b) 
Figure 4: Sample interactions of two massive particles. 

A different behavior is shown by the interactions of two massive particles, 
such as in Fig. |]a,b. The tree level result in our normalization is 

M tree = k 2 • (II.8) 

which again is dimensionless. Iterating this to form a loop gives us 

Mi oop ~ t£m\m\ ■ f dH — - 1 • 1 ■ 1 (II.9) 

J mi(l + p) m 2 {l+p) [l + q'Y [l + q) 1 

which by the same reasoning is 

Mioop ~ k ■ ~ n ■ ■ K 2 mi m 2 . (11.10) 

Here the expansion parameter appears as n 2 m 2 . An explicit calculation 
of this diagram later in this paper confirms that this is the correct result 
for the diagram by itself. This expansion parameter K 2 m 2 would cause the 
problem described in the introduction. 

Now let us turn to the general power counting result. Our goal is to obtain 
the power of q (with q being a typical external momentum) that a general 
diagram would yield. This will tell us what order in the energy expansion that 
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diagram will contribute to. The problematic class will be manifest by having 
diagrams with increasing number of loops which yield the same power of q, so 
that to calculate to this order in the energy expansion one would apparently 
needs to sum all the diagrams in this class. For a general result we need to 
allow for vertices not just from the lowest order gravitational action, but also 
from ones which contain more derivatives. Let us write this schematically as 

S g = fd 4 x^- \- \R + K 2 R 2 + 4R 3 + ..] (11.11) 

such that the coefficients of a gravitational Lagrangian with n derivatives 

will be -^5~- Note that Kq ~ cnergy . In a pure gravitational theory one 
would expect k ~ k, but there is no need to impose such a restriction here. 
Likewise the matter Lagrangian can involve extra derivatives on the light 
fields. We let the coefficients of the higher derivative terms involve a scale 
«o, i- e. 

S m = Jd 4 x^- (11.12) 

■ (dfldTQ - m 2 $ 2 ) + nfRd^d^ + 4 R2d ^ d ^ + ■■■ 

so that the coefficient of a Lagrangian with I derivatives on the gravitational 
field is m 2 Ko l . [Again, kq ~ and Tto can be kept distinct from k and k 
if desired.] 

Our procedure is to count powers of k, kq, 7to and m 2 in a general diagram. 
The remaining energy factor of the diagram, needed to give the proper over- 
all dimension, will come from factors of the external momenta. Consider a 
diagram with N% external matter legs and iVf external graviton legs, with 
a series of interactions between these particles. Correspondingly let Nf 1 and 
Nj be the number of internal matter and graviton propagators respectively. 
There are Ny vertices involving only gravitons, and N™ vertices which in- 
volve matter fields plus any number of gravitons, and a total of Nl loops. 
However, these vertices need to be categorized by the number of derivatives 
that are involved. For example, let iVy[n] be the number of graviton vertices 
which come from a Lagrangian with n derivatives. Clearly, N v = J2 n Ny [n] . 
Likewise the number of matter vertices with I derivatives on the light fields 
will be called N^[l] with iV{? = £,.MJ?[Z]. We illustrate this with a sample 
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Figure 5: Sample diagram with N 9 E = 1, N% = 4, Nf = 6, Nf = 4, 
N 9 V = 2, N{? = 6, N L = 3. 

diagram in Fig. ^. All matter lines propagate all the way through a diagram 
without terminating. 

With these definitions the coupling constants contribute the dimensionful 
factors 

(k 2 ) v -(k ) ■ (m) "(^) • (H.13) 

In addition, because each internal graviton line is formed using two vertices, 
the graviton fields will contribute a power of 

k) 1 E (11.14) 



from the normalization of the metric in (|II.2| ). On a matter line, there will 
be (v — 1) propagators if there are v vertices. Thus the number of matter 
propagators Nf 1 satisfies 

N™ = N^- -N%. (11.15) 

Since each propagator counts as a power of — , this contributes mass 
factors 

1 \ J V 2 "e 

(11.16) 

mj 
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These constitute all of the general dimensionful parameters except the 
external momenta and the loop momenta. When the loop integrals are reg- 
ularized dimensionally there will not be any powers of a regulator mass and 
the remaining dimensions after integration will be carried by the external 
momenta. Let us generically call these momenta q, and describe the power 
of the momenta by a factor q D . It is the dimension D which we are seeking 
in this exercise. 

Overall, this matrix element carries a dimension 

^-(Energy) E . (11.17) 

From our identification above, this is decomposed as 

/ \ 4_7V m — N 9 

A ~ (Energy) E E (11.18) 



^MKH f^\ 2N v f^\Ei l - N vW f^N 9 +N 9 E -2N, 



K ) ■ [mj ■ {Kq) 



1 \ "V 2 JV B 



'/??. 



There are however some relations among all the variables. For example, the 
total number of internal lines can be expressed in terms of the total number 
of vertices and the number of loops. The relation is 

N? + Nf = Nl + (Ny + Ny) — 1 (11.19) 

l n 

We can use this to eliminate Nf , using also NJ 1 = Ny — \N^, to find 

Nf = (Nf + Nf 1 ) - (n$ - -N%\ (11.20) 
= N L + lN% + Y,N v [n)-l. 



Plugging this into the general formula Eq. (|II.18|) , using J2n Ny [n] = N v and 
recalling the 
resulting in 



recalling that k, Kq, Kq all go as En ^ allow us to solve for the parameter D, 



Arm 

D = 2-^ + 2N L -N? + Y,{n-2)N$[n]+ y £l-Nv[t\- ( IL21 ) 

Z n l 
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This is our general power counting result for the momentum dependence 
of a general diagram. If we disregard the matter vertices, iVg = Ny[l] = 
Ny = 0, it is identical to Weinberg's theorem for chiral theories. The mo- 
mentum dimension of a diagram is higher if we increase the number of loops 
or if we use a gravitational Lagrangian with more than two derivatives. This 
shows that the power counting of loop diagrams in pure gravity involves the 
parameter K 2 q 2 (or K^q 2 if k ^ k.) 

In the presence of matter, the last term also behaves as expected: using 
the Lagrangian with extra derivatives on light fields (/ > 0) only increases the 
power of the momentum. However, the problem arises because of the minus 
sign in front of Nltf. Increasing the number of matter vertices in diagrams 
does not increase the order in the energy expansion. This cannot actually 
make the momentum power D decrease by increasing the number of matter 
vertices, since to add matter vertices to a given process we also have to change 
the number of loops. However, there are diagrams where one can increase 
the number of loops by one while increasing the number of matter vertices 
by two. Fig. ^ is one such example. This leaves D unchanged. Thus higher 
loop processes contribute at the same level to the energy expansion as tree 
processes. This gives a loop expansion of K 2 m 2 instead of n 2 q 2 . 

In summary, we have computed the momentum power or a given pro- 
cess ( |H.13| ), and found a class of dangerous diagrams where the addition of 
two matter vertices adds only one loop to the process, leading to no net 
increase in the momentum power D. 



Ill Analysis of a related model 

Because the couplings of general relativity are complicated by the tensor 
indices, it is easier to analyse a simpler model first. Although it is clear that 
the model is not identical with general relativity, it will nevertheless exhibit 
several interesting features which we will be able to generalize to the case of 
gravity in the next section. 

Consider a massless scalar field h coupled to one or more massive scalars 
$ with a trilinear coupling which carries the same strength as our counting 
rule given in ( |11.13| ). That is, we identify the Lagrangian 

£ = -df&dfQ + m 2 $ 2 (l + Kh) + ^hd^h. (III.l) 
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y%(m 2 ) ) ^4(^2) 

Figure 6: Tree level graph for heavy scalar scattering. 

The coupling nm 2 ^ 2 h enters into the power counting derivation in the 
same way as the lowest order graviton coupling and in this theory we obtain 
the same momentum power as in (|II.13|) with iVy[n] = for n > 2 and 
Ny [I] = for I > 0. The dangerous class of diagrams identified in the 
previous section also are equally problematic for this model. 

Let us verify the result of the counting theorem by considerations of the 
gravitational interaction of two heavy masses, as in Fig. [7]a,b. The single 
"graviton" exchange vertex, Fig. || has magnitude 

M = K 2 "^- (III.2) 

as expected. Note that in forming a nonrelativistic static potential one di- 
vides by 2mi • 2m,2 to account for our normalization of the states, obtaining 



cPq 1 2 mfm 2 



V{r) = iT^-T^—T^r-^^-e 



2 iq-r 



(27r) 3 2m! • 2m 2 q 
-J ■ " (III.3) 



Now consider the box diagram. Fig. [7|a: 

r d 4 k I 

M box = i ■ (nmj) 2 ■ (nml) 2 



(2tt) 4 [k 2 + ie][(q- k) 2 + ie] 

1 ' (ni.4) 



(pi + k) 2 — m 2 (p 2 — k) 2 — ttt-2 
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Pi{mi] 



k» 




Jq^-k" 



(a) 



(b) 



Figure 7: The (a) box and (b) crossed box graphs which have the wrong 
naive power counting behavior. 



We combine the graviton propagators using the usual identity 
11 r l dx 



k 2 (k + qf 



o [k 2 — 2xk ■ q + x 2 q 2 } 
1 dx 



(k — xq) 2 — Q 2 



(III.5) 



with Q = —x(l — x)q 2 , and do likewise for the matter propagator. The 
diagram is ultraviolet finite but has an infrared divergence which we regulate 
dimensionally. Integrating in d dimensions we obtain 



M box 
with 



T 4 - - 

[Km\f ■ {Kmlf ■ (j, 4 *-* 1 2 



(4tt) 



dx f dy jT dz- V ^_d 



o 



(III.6) 
(III.7) 



Mi 



z m 1 + (1 — z) m 2 — 2z(l — z)p\ ■ p 2 - 



For small — , the y integration can be done near d = 4 dimensions, 



o L> 4 ~2 



2M 2 



Q 



6~d 



1 + {d - 4) 



1 Q s 



m 



(III.8) 



The infrared divergence as d —>■ 4 is in the x integration. In the z integration 
there is a pole when the intermediate matter state is on shell and this is made 
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well defined by the usual ie factors. However, both of these integrals can be 
done straightforwardly. If we define p\ ■ P2 = m\m<i + w, we end up with 



M 



box 



Q 



A-d 



16vr 2 



In 



w 



+ m- 



3mim.2 p(mi + 1712) 
+ constant + O 



(III.9) 




with p —\ p I in the center of mass. If we defer comment on the imaginary 
part of this amplitude to below, we see that this does obey the expectation 
of the power counting theorem, with a correction of order /? 2 m 1 m 2 compared 
to the tree level amplitude. 

However, we must also consider the crossed box diagram in Fig. [7]b. This 
can be handled in a similar fashion, and is slightly easier because it does not 
have an imaginary part. Defining p\ ■ p 4 = m\m<i + w', we obtain 



M 



crossed 



K 2 rn\rri2 K 2 mim 2 



A-d 




+1 



w 



2>m\m2 



(111.10) 



^-r ] + constant + O ( 
2 ' \M , 



rn 



with the same constant as in (|III.9|) . We see that the most dangerous terms 



cancel between the diagrams. Using w 



■w 



we get the final result 



M 



total 



Mtree + 

2 2 2 



■^^■box ■A" 1 ! crossed 
1 



;m.n) 



1 + 



16% 
2 



A-d 



1 2 2 K m\V(l2 

-K q +Z7T— ^ 

p[mi + 1712] 



constant 




In the real part of the amplitude, the expansion parameter has become 
K 2 q 2 , which is well behaved. The imaginary part of the amplitude is just a 
phase and does not contribute to observables at this order when the matrix 
element is squared. It is simply the analogue of the well known "Coulomb 
phase" and is generated from the rescattering of the on shell intermediate 
state of matter particles. Like the Coulomb phase, it has been shown by 
Weinberg pj to exponentiate to all orders in general relativity (the proof 
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extends to this simpler theory as well), so that this term does not cause any 
trouble. The same paper by Weinberg also proves that infrared divergencies 
cancel in general relativity and by extension in this theory, by the considera- 
tions of virtual corrections of these diagrams plus the bremsstrahlung radia- 
tion of real particles. These can be regulated dimensionally also ||, and will 
yield as finite effects residual corrections of order K 2 q 2 and K 2 q 2 Inq 2 . We are 
not here directly interested in the exact answer; for us the important result 
was the cancellation of n 2 m 2 effects in ( |IH.11|) . 



V 



p+q 



p+q 




V 



p+q 




(a) (b) (c) 

Figure 8: The definition of the 'Bose-symmetrized' propagators. 

We can show that this cancellation is not peculiar to the box diagram, 
but rather is a general feature of this theory. This can be seen by considering 
the basic unit of a single line with two interactions with off-shell gravitons, 
as in Fig. [8[ Any time a given ordering is possible, the crossed order is also 
possible, cf. Fig. |8|a and [8|b. If we add these two diagrams and allow the 
external legs to be on-shell, we find that the sum of propagators behaves as 
— s whereas each individual propagator was of order — : 



V 



Km 



[Km 



[Km 



1 



+ 



(p + k) 2 — m 2 (p' — k) 2 — m 2 



2\2 



+ 



2\2 



2p-k + k 2 -2p' -k + k 2 

2{p' - p) ■ k - 2k 2 
[2p-k + k 2 } ■ [2p' -k-k 2 ] 



Km 



(111.12) 



However, since p' — p = q, there is no factor of the large mass m in the 
numerator, 

q 2 -k 2 -(k- q) 2 



2\2 



V = {Km 2 ) 



[2p-k + k 2 } ■ \2p' -k-k 2 } 



(111.13) 
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Because there are two factors of p ■ k ~ m in the denominator, this double 
vertex counts as 

V ~ ^ ~ K 2 m 2 (111.14) 

rather than the V ~ n 2 m 3 that the naive counting would imply. For Fig. [8]a or 
individually would give the extra two factors of — which converts the unde- 
sirable expression to n 2 q 2 , thereby explaining the result found above. 

The only exceptions to this power counting occurs for what can be termed 
"exceptional momenta." This refers to momenta where the propagator is 
not of order — , and can occur when the intermediate line goes on shell, 
(p + k) 2 — m 2 = = 2p ■ k + k 2 so that p ■ k ~ k 2 . In this case we do 
not gain a power the power of — from the propagator. This is exactly what 
was found above in the explicit calculation of the box diagram. The on-shell 
intermediate states generate the imaginary part of the diagram which has a 
different dependence on the masses than does the remainder of the diagram. 
This leads to the Coulomb phase in the box diagram. 

p p+q 



T>Zi+§ T|/ 2 +f T|ijf§ 

Figure 9: The two- loop 'Bose-symmetrized' propagator, defined for J2\ h = 
0. 

We have been able to extend the demonstration of cancellations to three 
and four vertices (using computer algebra). The difficulty here appears be- 
cause only the external lines are on shell. The desired cancellation does not 
occur for just the symmetrized sum of any two of the permutations, e. g. when 
we permute l\ and I2 in Fig. [| but requires the sum of all six permutations. 

We take the 'Bose-symmetrized propagator' denoted by a double line 
in Fig. ||: 

D^(p,h,l 2 ,l 3 ) = ^. £ 



>f) 2 -m 2 



3' p&miP + h + D 2 -m 2 (p + h + l 2 + 2\ 

(111.15) 

where the summation goes for all permutations of the k, subject to J2i h = 0. 



14 



Naive power counting would say ~ O (j^p) an d for the total dis- 
appearance of the kitl expansion parameter we would need ~ O 
because the expected behavior of the amplitude is 

in power counting: 



' I 1 + Mone-loop + (^two-loop + ■ ■ } ( IIL16 ) 



and in the effective theory 



q 2 



{l + Mone-loop + (^)two-loop + •••}• ( IIL17 ) 



A tedious but straightforward algebra shows indeed that, disregarding 
special points in the integration region where some of the denominators are 
less than 0(rrii), cancellations occur with the result 

D^(p,h,l 2 ,l 3 ) = 0^y (111.18) 

Although we do not have an inductive proof valid for all orders we did 
the same calculation for D$(p, h,l2,h, U), figuring in ladder- type three- loop 
diagrams, which indeed showed a similar behavior 

D^\pMMMM) = 0[^j (111.19) 

in accordance with the requirement of the validity of the energy expansion. 
This calculation was done with a Mathematica program and attempts to 
calculate D$ did not succeed because of the prohibitively large amount of 
CPU time required. We have unfortunately not found an inductive proof 
that lets us extend these results to all orders. These calculations show that 
in the toy model, which is free of gauge complications, the power counting 
theorem works correspondingly to Weinberg's power counting theorem in 
chiral effective QCD due to unexpected cancellations. 



IV Lessons for general relativity 

Many of the features uncovered in the last section are also applicable for gen- 
eral relativity. Although general relativity contains additional multigraviton 
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vertices, the matter vertex considered in the previous section is the most 
dangerous. We have not been able to reformulate the general power count- 
ing formula (in harmonic gauge) such that the cancellation is explicit in all 
diagrams. However, we can display how the desired behavior is restored in 
examples which are physically relevant. 




(a) (b) 
Figure 10: Graviton vertex and one of the loop corrections. 

Before we turn to the box diagram once again, let us consider the vertex 



correction, Fig. [10]. Of the several contributions to the vertex at one loop 
only Fig. [I0[ b and the self-energy diagrams are of the dangerous category, 
with two extra matter vertices and one loop. 

Because the vertex coupling is the energy momentum tensor, and the 
energy momentum tensor is conserved, there is a nonrenormalization theorem 
for the matrix element at q 2 = 0. The general form for the vertex, consistent 
with the conservation d^T^ v = is 

< p' | T M „ | p >= Fx(q 2 )- [pftpl + p v p'^ - 7^] +F 2 (q 2 )-[q fl q u - ij^q 2 ] (IV. 1) 

and at tree level F\ = 1, F 2 = 0. 

The dangerous diagrams naively give a correction to F% of order K 2 m 2 , 
and individually will do so. [Because powers of q 2 multiply F 2 it automati- 
cally is not a problem - a one loop contribution to F 2 (q 2 ) of order K 2 m 2 is 
allowable in a well behaved energy expansion.] However, since T^ v measures 
the physical energy and momentum we have the constraint -Fi(O) = 1. Thus 
all contributions to Fi(q 2 ) which are independent of q 2 , in particular all cor- 
rections of order K 2 m 2 , must cancel when expressed in terms of the physical 
mass and momenta. This occurs by the cancellation of the vertex pT0] b with 
the renormalization due to the self energy. This is entirely analogous to the 
nonrenormalization of the charge form factor in QED at q 2 = 0. 

One can also repeat the exercise to show that the sum of the two di- 
agrams in Fig. [7]a,b, but with two gravitons instead of scalars, behaves as 
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K 2 m 2 instead of the K 2 m 3 behavior as given by the naive power counting. 
The gravitational vertex is more complicated than the scalar one, for exam- 
ple involving np^p + k) u at a given vertex instead of urn 2 . However, in the 
counting of powers of mass p^k v is already one factor fewer power of the mass 
than is p^p u . Therefore in showing that the K, 2 m 3 behavior is not present in 
the sum of diagrams, we need only consider the np^Pu portion of the vertex, 
which is common to all vertices and which will not upset the cancellation of 
the two diagrams. The proof then goes through exactly as in the last section. 

Finally, we seek to demonstrate the cancellation in the box plus crossed 
box diagrams. 

In order to show that the situation, including the necessary cancellations, 
is analogous to the one in the toy model, we calculate the contribution from 
the box and crossed box graphs on Fig. [F]a,b and we will see that the can- 
cellation of the > O Mj^ terms does occur in the real part of the amplitude 
and that in the nonrelativistic limit, i. e. for small values of spatial momenta 
p <C y/s, the toy model gives the correct coefficients of the leading terms. 
This could be expected as we will explain later. 
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i 
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P P' 

| => -|« • (p, pi) 

Figure 11: Propagators and the relevant vertex in harmonic gauge. 
In harmonic gauge the graviton propagator is 

a 3 1 nB 



■ P map with P map = ~ , (IV-2) 



the heavy mass propagators are 



p2 _ m 2 



(IV.3) 
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and the graviton-scalar-scalar vertex is 

-in ■ V™ (p, p') with (p, p') = p^p' v) - -(p-p'-m 2 ), 

(IV.4) 

see also Fig. [TT| Using these rules and the usual Feynman parameterization 
in the loop integral we arrive at 
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(IV.5) 



2-r i - 



ox 



+ 



+ r|-- 



(c)box 



Tp/j,ua/3 
''"■I ' r (c)box 



(c)box 



where the IR singularities (at e —>■ 0) come from the first term and the UV 
singularities come from the last term. The functions F are determined as 



box 



cbox 



Vpfapi + k) ■ pr<# . V$(p 2 ,p 2 - k) 

V^(pi +k, Pl + q)- P x ^ 6 ■ V^(p 2 -k,p 2 - q) 
Vp(pi, P i + k) ■ pr<# ■ V^(p 2 - q + k,p 2 -q) 
V x 7(Pi +k, Pl + q)- P xp > KS ■ V^(p 2 ,p 2 -q + k) 



(IV.6) 
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(IV.7) 
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u = 2{m 1 + m 2 ) — s — q = Uq — q . 



(IV.8) 
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(a) (b) 

Figure 12: (a) The complex path of z integration, (b) What this path 
becomes in the complex y plane. Singularities are shown in red, the path of 
integration in green. 



The quantity [ }f ox is always positive, corresponding to the fact that the 
crossed box diagram has no imaginary part (which would emerge from the 
point where M?^ box — > 0). The corresponding [ ] b ° x does have zeroes which 



can be avoided by a complex continuation in z as shown in Fig. 12 



After some a tedious algebra we managed to separate out the IR divergent 
part which only stems form the x = H,y = corners of the integration 
region. The coefficients were extracted with a REDUCE program and cast 
in the form 



\k=l 



iM {c)box = - 8 ( 4?r p J Q dz J Q -cos (3 sin p ■ \^J2hk(l,z)ipk(b) J (IV.9) 
with 

cos (3 + sin /3 



7= cos/3 -sM 



and 



and 



[ - cos 2 P+[) z -sin 2 fi-ie, (IV.ll) 



Mb) = \, Mb) = ^, Mb) = In 6, Mb) = (IV.12) 
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and the functions hk(j,z) turn out to be polynomials in q 2 ; they contain 
about 932 terms that we do not display here. Their algebraic structure can 
be understood though and we find that the only nonpolynomial dependences 
in any of the variables come from seven special functions of 7 and their 
singularities are avoided by the above prescription to use a complex path of 
integration for z. 

The calculation at this point contains no approximation at all. However, 
the integrals cannot be done by quadrature. For our purposes it is only nec- 
essary to separate the terms which are divergent in the limit q — * 0, and due 
to the fact that the integration contours defined above avoid the singularities 
we can bring the q — > limit into the integrands. For lack of space we do 
not describe all the details of the thorough analysis of the divergent terms. 
What we do is we first observe that cutting the f3 integration in two parts, 
jp( s , mi ,m 2 ) _|_ ^ e j n t e g ranc i j n the latter is uniformly bounded and so 
provides no divergent terms. Using a Taylor- formula for the former part, the 
remainder is also seen to be uniformly bounded (this was proven for each k 
using Weierstrass's theorem) and what remains is 




(IV.13) 





([ ]z) q ^ and h (0) (x,z) = /i(x->tan/3, z) q ^ Q . 



lim^o 4 2 4 and the prime always refers to the derivative 



20 



in x and we used the notation 

4>{z) = { 



1 + 4z(l — z) for cbox 
1 for box. 



Here, h 4 =8 [s 2 — 2s{m\ + m 2 ,) + (mf + mf)] is independent of /3, z. 
This term generates the leading divergence, and as we will see this term 
is the only one present in the toy model of Sec. |TTT| . 

The most divergent terms are multiplied by the integrals Jq -%y- A sub- 

[ ] z 

stitution of variables z — >y with - + - = 2, has the interesting property 



]box 



r -i (O,cbox) 

(2y-l) 5 



and tracing the path of the complex y integration as shown in Fig. 
find a residue term in 
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(IV.14) 



i.e. the sum of the two integrals figuring in (|IV.13|) is a pure residue. We 
find that the leading divergence is IR and UV finite: 

.1 



lng 



K 



167T q 2 py/s 



2s(m 1 + m 2 ) + (m 1 + m 2 ) 



(IV.15) 



is pure imaginary and as such can be included in an unobservable phase. 

The next-to-leading, \ divergence contains the same factor and also re- 
sults in a pure imaginary contribution 



16n q 2 py/s 



s 2 — 2s{m\ + m\) + (m 1 + rri. 




(IV.16) 



where the e — > singularity is of IR origin. 

The ambiguity in ( IV.15 ) in defining In q = In a + ln/i is a reflection of 
the ambiguity in the IR regulator to remove IR divergencies. Total removal 
is achieved only after taking into account soft graviton processes as in J4j and 
different choices of fi show up as different constants in ( [1V.16| ): 
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(IV.17) 
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The O (~j divergencies are much harder to calculate. We found, by a 
similar calculation, that the O part has no imaginary part, contains no 
IR or UV singular terms. Its coefficient is a complicated function of s, m\ 
and r»2 which we do not display for lack of space but whose p — > limit is 



4 2 2 

k m\m,2 
128g 



mi + m 2 + O ( — 



rn 



(IV.18) 



This term, which still contains one power of nm too much, corresponds 
to a contribution to the classical potential 



V(r) — -G^ 



G(mi + m 2 ) _^ 
4r 



(IV.19) 



and forms part of the classical general relativity corrections which are en- 
coded in the one-loop diagrams of Fig. |7] when the harmonic gauge is used. 
We note that no similar real term arises in the toy model of Eq. fllll.lip . 

The p — > limit of the coefficients of q — > divergent terms turn out to 
be smooth except for the imaginary part which has a K 4 m term in ( |1V.16| ). 
The limit we use in this calculation, 

— > - -> (IV.20) 

m p 

exactly corresponds to the requirement that no bound states are formed 

p 2 p 2 
Gm\m2 ■ q ~ h 



2m\ 2m 2 

i. e. (IV.21) 

p 2 

2Gmii7i2 > oo. 

q ■ m red 

We finally comment on why out toy model gave the correct leading di- 
vergencies. At leading order in — , the heavy scalar vertex becomes 

-7^ ■ W ' Pu+ Vv ■ p'p- Vvv ■ (P-P ' - m2 )} — ► -iKm 2 v^v u . (IV.22) 

This vertex is attached to a graviton propagator and is multiplied by P£g 
and so becomes ^ 

-iKvn 2 ■ (v a vp - yVatp) ' (IV.23) 
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so when we multiply this by the heavy scalar on the other end of the graviton 
propagator, the full theory, as compared to the toy model, provides and extra 
factor 



(vi • v 2 ) 



(2pi ■ p 2 ) 2 — 2m\m\ _ [s 2 — 2s{m\ + m|) + (m\ + mf)] 



2 Am{m2 Amlm^ 

(IV.24) 

Comparing now (PV15D with the appropriate, O (^) part of ( PTTHl) wc sec 
that the square of ( |1V.24| ) shows up, corresponding to the presence of two 
graviton lines. A similar comparison between flTV.lSfl and flM.ll|) shows that 
in the toy model the ^ divergence is absent. It is nonzero in ( [L V . 1 8| ) only due 
to subleading contributions to (|IV.22|) and (|IV.23|). 



V Power counting in a physical gauge 

In covariant gauges, both classical and quantum effects are included in the 
same Feynman diagram. The simplest example is the one graviton exchange 
diagram, which includes the classical Newton potential, but loop diagrams 
also exhibit this property 0. However, in physical gauges, classical and quan- 
tum effects are separated. The physical quantum degrees of freedom are trans- 
verse and traceless 0], corresponding to massless spin-two quanta. In a mul- 
tipole expansion the monopole term, which generates the Newton potential, 
is classical, while the spin-two degrees of freedom couple to the quadrupole 
term. This suggests that the dominant urn 2 coupling that caused us do much 
trouble in Sec. |TJ is to be associated with classical physics while the quantum 
degrees of freedom have a milder behavior. We will show that this is the 
case for the interaction of two heavy masses. This allows a quantum power 
counting which is well behaved. 

Covariant gauges, in particular the harmonic gauge treated covariantly, 
are preferred for practical calculations 0. When combined with the back- 
ground field technique, they can explicitly retain the invariances of general 
relativity. In contrast the construction of the physical gauge quantum the- 
ory severely disturbs the underlying symmetry of general coordinate invari- 
ance |5|. One picks a preferred frame for the quantization, and the split of 
quantum and classical degrees of freedom depends on that frame. Never- 
theless a physical gauge is attractive conceptually because only the physical 
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radiation degrees of freedom are quantized. This is analogous to the Coulomb 
gauge quantization of QED. 

Although the general metric tensor has ten real components, the 
radiation field has only two independent degrees of freedom, corresponding 
to helicity ±2 ||. For a wave propagating in the z direction, a harmonic 
gauge constraint plus residual gauge freedom can be used to reduce the ten 
original components of a polarization tensor e^ u such that only en = — e 2 2 
and 612 = 621 are nonvanishing (other choices are also possible.) The helicity 
states with A = ±2 can be found from these via 

e± = en =F ^12 = -e 2 2 =F ^12- (V.l) 

When quantized, these become the graviton degrees of freedom. 

Let us consider the interaction of two heavy masses, which have a small 
relative momentum. The momentum of one of the masses can be written as 

P M = mv^ + p^ (V.2) 
v„ = (1,0,0,0V 




(a) (b) 
Figure 13: Suppressed matter-gravity interaction vertices. 

Since 

%=V P ,- %, (V.3) 

we will treat the residual p as of order q in the energy expansion. If we 
quantize in the frame where defines the timelike direction, then the phys- 
ical graviton degrees of freedom will only couple to the spacelike traceless 
components of the energy- momentum tensor T^-. Since the large mass does 
not contribute to these, we know 

Tij^P^+p'jPi (V.4) 
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and we count as order q 2 , whereas Too is of order m 2 . The same suppres- 



sion is present for the two graviton interaction of Fig. |13|a, and the coupling 
of one classical Newtonian field and one physical quantum field as shown in 
Fig. Ob. From the general vertex 



K 



- - • {j]rj\ ■ Ifac0 + Vpcr ' Ir,\, a p) ' p'V (V.5) 



we have 



Too.oo ~ K 2 m 2 



T 00,ij 



2 2 

k q 

2 2 

k q . 



(V.6) 



What we see from these couplings is that the physical quantum degrees of 
freedom have a reduced matter coupling with two fewer powers of m com- 
pared to the harmonic gauge counting rules. 




(a) 



(b) 



(c) 




(d) W 
Figure 14: Interaction of two masses in a physical gauge. 
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Various other diagrams describing the interaction of two masses in phys- 



ical gauges are shown in Fig. |1J. The dashed line represents the classical 
Newtonian potential while the wavy line describes the quantum degrees of 
freedom. In the normalization used throughout this work, the classical in- 

2 4 

teraction is again of order ^Mjr-. Corrections due to vacuum polarization 
in Fig. |H]p,c have the same powers of m but two extra powers of k 2 , leading 
to a result of order 

2 4 

Km 



, K 2 q 2 -(a + b\ogq 2 ) (V.7) 
T 

so that these diagrams have the expected expansion parameter n 2 q 2 as de- 
sired. An important feature of the physical gauge is that diagrams with 
extra matter couplings are suppressed. For example, the mixed box diagram 
of Fig. [I4|d has two factors of — from the propagators, but no compensating 



factors of the mass in the vertex coupling or the physical gravitons. This 
means that this diagram is suppressed by ^ compared to the vacuum po- 
larization corrections. All of the diagrams with graviton-matter vertices are 
suppressed by powers of This leaves the vacuum polarization diagram as 
the leading quantum correction is this gauge, with a well behaved expansion 
parameter. 

The above results may be seen more easily in different normalization for 
the fields, typical for the nonrelativistic limit. Here, we divide all matter 
vertices by a factor of 2m, so that T Q0 = m for a particle at rest. In this 
normalization (used also in Heavy Quark Effective Theory |J) propagators 
have no factor of — . However, the coupling of the transverse traceless degrees 
of freedom to matter fields are proportional to 

Ta ~ 2m (V.8) 

and vanish as m — > 0. Then all diagrams with matter coupling of graviton 
simply drop out, leaving only the vacuum polarization diagram in this gauge. 

To close this section, we note that the use of a physical gauge seems to 
be required if we are to be able to introduce the idea of a purely classical 
source. In the harmonic gauge, the inclusion of the matter propagator was 
required to properly identify the classical and quantum corrections. Both 
vertex and vacuum polarization diagrams are important. One could not at 
the start of the calculation take the mass to infinity and treat this resulting 
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field in the classical limit. The reason is that the vertex coupling strength 
also goes to infinity in this limit. However, in the physical gauge, the dia- 
grams with matter propagators and couplings to transverse quantum fields 
are unimportant in the limit m — > oo. By taking this limit in this gauge, 
one obtains a classical source, with an interaction which receives quantum 
corrections. 

VI Summary 

We were motivated for this study by the observation of a class of diagrams 
which in harmonic gauge would apparently upset the utility of the energy 
expansion, and indeed also spoil the classical limit of the theory. Part of the 
problem is due to the fact that the graviton propagator in harmonic gauge 
includes both classical and quantum effects. By consideration of several 
diagrams we were able to demonstrate the nature of the cancellations in 
harmonic gauge which removed the undesirable expansion parameter and 
led to a well defined energy expansion. The logic for this behavior is clearer 
in a physical gauge, analogous to Coulomb gauge in QED, even if explicit 
calculations are much more painful in such a gauge. The physical transverse 
traceless quantum degrees of freedom only couple with a reduced strength in 
problems with nearly static matter fields. In the limit that m becomes very 
large, the effect of matter couplings become negligible, and the modification 
to the gravitational self-interactions (i. e. vacuum polarization) become the 
most important quantum corrections. These self-interactions of the sector 
satisfy the Weinberg power counting theorem without any problems. This 
indicates that the quantum energy expansion is well behaved in physical 
gauges, and hence by extension in all gauges as long as one is calculating 
gauge invariant quantities. 
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